Time-dependence of correlation functions following a quantum quench 
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We show that the time-dependence of correlation functions in an extended quantum system in d 
dimensions, which is prepared in the ground state of some hamiltonian and then evolves without 
dissipation according to some other hamiltonian, may be extracted using methods of boundary 
critical phenomena in d + 1 dimensions. For d = 1 particularly powerful results are available using 
conformal field theory. These are checked against those available from solvable models. They may 
be explained in terms of a picture, valid more generally, whereby quasiparticles, entangled over 
regions of the order of the correlation length in the initial state, then propagate classically through 
the system. 
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Suppose that an extended quantum system in d dimen- 
sions (for example a quantum spin system) , is prepared at 
time t = in a pure state \ipo) which is the ground state 
of some hamiltonian Ho (or, more generally, in a ther- 
mal state at a temperature less than the gap mo to the 
first excited state.) For times t > the system evolves 
unitarily according to the dynamics given by a different 
hamiltonian H, which may be related to Hq by varying 
a parameter such as an external field. This variation, or 
quench, is supposed to be carried out over a time scale 
much less than m^ 1 . How do the correlation functions, 
expectation values of products of local observables, then 
evolve? The answer to this question would appear to de- 
pend in detail on the system under consideration. It was 
first addressed in the context of the quantum Ising-XY 
model in Refs. (see also 0). Until recently it was, 
however, largely an academic question, because the time 
scales over which most condensed matter systems can 
evolve coherently without coupling to the local environ- 
ment are far too short, and the effects of dissipation and 
noise are inescapable. However, with the development of 
experimental tools for studying the behavior of optical 
lattices of ultracold atoms, and quantum phase transi- 
tions in these systems , there has been renewed interest 
in this theoretical problem (see, for example, 0,0-) 

In this Letter we study such problems in general and 
argue that, if H is at or close to a quantum critical point 
(while H is not), there is a large degree of universality in 
the behavior at sufficiently large distances and late times, 
despite the fact that correlations typically fall off expo- 
nentially, rather than the power laws characteristic of the 
ground state near a quantum critical point. Our argu- 
ments are based on the path integral approach and the 
well-known mapping of the quantum problem to a classi- 
cal one in d+1 dimensions. The initial state plays the role 



of a boundary condition, and we are able to then use the 
renormalisation group (RG) theory of boundary critical 
behavior (see, e.g., |fj). From this point of view, particu- 
larly powerful analytic results are available for d = 1 and 
when the quantum critical point has dynamic exponent 
2=1 (or, equivalently, a linear quasiparticle dispersion 
relation u> = v\k\) because then the 1 + 1-dimensional 
problem is described asymptotically by a boundary con- 
formal field theory (BCFT)0,@. Some of these methods 
have recently been applied |fj to studying the time evolu- 
tion of the entanglement entropy, but, as we shall argue, 
they are more generally applicable. Further details of 
these calculations will appear elsewhere [in| . 

The results we find from CFT suggest a rather sim- 
ple picture which is, however, more generally applicable: 
the state |^>o)i which has an (extensively) high energy 
compared with that of the ground state of H, acts as 
a source for quasiparticle excitations. Those quasipar- 
ticles originating from closely separated points (roughly 
within the correlation length £o of the ground state of Ho) 
are quantum entangled. However, once they are emit- 
ted they behave semi-classically, travelling at speed v. 
They have two distinct effects. Firstly, incoherent quasi- 
particles arriving a given point r from well-separated 
sources cause relaxation of (most) local observables at 
r towards their ground state expectation values. (An ex- 
ception is the local energy density which of course is con- 
served.) In the CFT case this relaxation is exponential 
~ exp(— Trxvt/2To). Here To ~ t 71 ^ 1 is non-universal, but 
x, the bulk scaling dimension of the particular observ- 
able, is related to the critical exponents of the quantum 
phase transition, and is universal. Hence ratios of decay 
constants for different observables should be universal. 
Secondly, entangled quasiparticles arriving at the same 
time t at points with separation |r| 3> £o induce (quan- 



turn) correlations between local observables. In the case 
where they travel at a unique speed v, therefore, there 
is a sharp 'light cone' effect: the connected correlations 
do not change significantly from their initial values until 
time t ~ |r|/2u. In the CFT case this light-cone effect 
is rounded off in a (calculable) manner over the region 
t— |r|/2i> ~ To, since quasiparticles remain entangled over 
this distance scale. After this they rapidly saturate to 
time-independent values. For large separations (but still 
■C 2vt), these decay exponentially ~ exp(— 7ra|r|/2vTo). 
Thus, while the generic one-point functions relax to their 
ground-state values, the correlation functions do not, be- 
cause, at quantum criticality, these would have a power 
law dependence. Of course, this is to be expected since 
the mean energy is much higher than that of the ground 
state, and it does not relax. 

These results rely on the technical assumption that the 
leading asymptotic behavior given by CFT, which applies 
to the euclidean region (large imaginary times), may sim- 
ply be analytically continued to find the behavior at large 
real time. While such procedures have been shown to 
give the correct behavior for the time-dependent corre- 
lations in equilibrium, it is important to check them in 
specific solvable cases. We have done this for the case of 
a lattice free boson (coupled harmonic oscillators), and 
for the Ising-XY chain, which can be transformed into a 
free fermion problem. These models, which are exactly 
solvable on the lattice, confirm our general results and 
also show how the semi-classical picture discussed above 
is modified for a more general quasiparticle dispersion re- 
lation u> = fifc , taking into account both the effects of the 
lattice and of a finite gap. The results are consistent with 
this picture as long as the quasiparticles are assumed to 
propagate at the group velocity fl' k appropriate to their 
wave number k. In this case the light-cone effect first 
occurs at time t ~ |r|/2?; m , where v m is the maximum 
group velocity. (If this occurs at a non-zero wave num- 
ber, it gives rise to spatial oscillations in the correlation 
function.) However, because there are also quasiparti- 
cles moving at speeds less than v mi the approach to the 
asymptotic behavior at late times is less abrupt. In fact, 
for a lattice dispersion relation where f2' fc vanishes at the 
zone boundary, the approach to the limit is slow, as an 
inverse power of t. A similar result applies to the 1-point 
functions. This is consistent with the exact results found 
inQ. 

We first discuss some general features of the path in- 
tegral approach and its relation to boundary critical be- 
havior. The expectation value of any product of local 
operators is given by 

(0(t,{Ti})) = Z- 1 (^|e« r *- eH 0({r i })e- <fft - eH |V>o) 

(1) 

where we have included damping factors e eH in such a 
way as to make the path integral representation of the 
expectation value absolutely convergent. The normal- 



ization factor is Z = (i/jo\e 2eH \4>o). Eq. JIJ may be 
represented by a path integral in imaginary time r 

i J [#(r,r)]O({r i },O)e- s ^^o|0(r,T 2 ))( ( /)(r,r 1 )|^o> 

(2) 

over a complete set of fields 0(r, r) (or, in a spin sys- 
tem, a coherent state representation), with S — J T ^ Ldr, 
analytically continued to t\ = — e — it and T2 = e — it. 
Here L is the (euclidean) lagrangian corresponding to the 
dynamics of H . We will consider the equivalent slab ge- 
ometry between r = and t = 2e, with O inserted at 
t = e + it. 

Eq. J2J has the form of the equilibrium expectation 
value in a d + 1-dimensional slab geometry with partic- 
ular boundary conditions. We wish to study this in the 
limit when t and the separations |r^ — Yj | are much larger 
than the microscopic length and time scales, when RG 
theory can be applied. If H is at or close to a quantum 
critical point, the bulk properties of the critical theory 
are described by a bulk RG fixed point (or some rele- 
vant perturbation thereof). In that case, the boundary 
conditions flow to one of a number of possible bound- 
ary fixed points 0- Thus, for the purpose of extracting 
the asymptotic behavior, we may replace |?/>o) by the ap- 
propriate RG- invariant boundary state |^q) to which it 
flows. The difference may be taken into account, to lead- 
ing order, by assuming that the RG-invariant boundary 
conditions are not imposed at r = and t — It but at 
t = —To and t = 2e + tq. In the language of boundary 
critical behavior, To is called the extrapolation length |6|]. 
It characterizes the RG distance of the actual boundary 
state from the RG-invariant one. It is always necessary 
because scale-invariant boundary states are not in fact 
normalizablc 8]. It is expected to be of the order of the 
typical time-scale of the dynamics near the ground state 
of i/o, that is the inverse gap m^ 1 . The effect of in- 
troducing To is simply to replace e by e + To. The limit 
e — ► 0+ can now safely be taken, so the width of the slab 
is then taken to be 2ro- 

d = 1 and CFT. The above was completely general, but 
we now consider the case when H is at a quantum crit- 
ical point whose long-distance behavior is given by a 
1+1-dimensional CFT, with dispersion relation uj = v\k\. 
We set v = 1 in the following. RG-invariant bound- 
ary conditions then correspond to conformally invariant 
boundary states. Under these circumstances the corre- 
lation functions of local operators in the slab geometry 
(whose points are labelled by a complex number w with 
< Imw < 2to) are related to those in a half-space 
Imz > by the conformal mapping w — (2tq / tt) log z . 
In the case where O is a product of local primary scalar 
operators <&i(wi) (which covers most of the physically rel- 
evant cases, for important exceptions see below) we have 



i i i 

where Xi is the bulk scaling dimension of Note that 
the expectation values of the <&i in the ground state of 
H are supposed to have been subtracted off. We now 
discuss some special cases of 

One-point junctions. In this case scale invariance implies 
that ($(z))uhp = A*(2Imz)~ x , where Af is an ampli- 
tude depending on the selected boundary condition, (ft 
may vanish if $ corresponds to an operator whose expec- 
tation value in \iPq) vanishes.) This gives 



($H) 



strip 



= Al 



i 



4r sin(7rr/(2r )) 



Continuing to r = tq + it we then find 



<*(*)} = Af 



I 



4r cosh(7rf/(2T )) 



(4) 



(5) 



which exhibits exponential decay for t ^> tq, with a decay 
time related to the critical exponent X, (IS described in 
the introduction. 

An important exception to this law is the local en- 
ergy density (or any piece thereof). This corresponds to 
the tt component of the energy-momentum tensor T^. 
In CFT this is not a primary operator. Indeed, if it 
is normalised so that (T^uhp = 0, in the strip JJjj 
(T tt (r, r)) = 7tc/24(2to) 2 (where c is the central charge 
of the CFT) so that it does not decay in time. Of course 
this is to be expected since the dynamics conserves en- 
ergy. A similar feature is expected to hold for other local 
densities corresponding to globally conserved quantities 
which commute with H , for example the total spin. 
Two-point functions. The 2-p_oint function in the half- 
plane takes the general form 



UHP — 



Z Y1 Z 2\ 



z 12^12 z ll 2; 22 



F(V) , (6) 



where z, 



z\ is the image of z, in the real 



axis, and r\ = ^ii^s/^iS^al- The function F is universal 
but depends in detail on the particular BCFT. Using 10, 
analytically continuing, and assuming r, t 3> tq we find 



($(Q,t)*(r,t)) = (V2r ) 5 



where now 



V ■ 



r/2r 



= 7rt/T 



e nr/2T _|_ g7rt/r 



(7) 



(8) 



Thus for r - 2t > r , 77 - e ^(t-r/2)/r < 1; while for 
2t — r To, 77 ~ 1. In both of these limits the behavior 
of _F is determined by the short-distance expansion. The 



first limit, when t < r/2, corresponds in the half-plane 
to both points being close to the boundary, when we can 
use the bulk-boundary expansion to argue that F(rf) ~ 
(Af) 2 r) Xb , where Xb is the boundary scaling dimension of 
the leading boundary operator to which $ couples. This 
gives, for t < r/2, 

($(r,t)*(0,i)) - (Af) 2 e-" x V To x e -***W-*V-* . (9) 

Note that if ($) 7^ 0, x b = and © is just ($) 2 . In 
that case the leading behavior of the connected 2-point 
function is given by subleading terms either in F or in 
the bulk-boundary short-distance expansion. The oppo- 
site limit, when t > r/2, corresponds to the two points 
both being far from the boundary, in which case the bulk 
behaviour dominates and F — > 1. This gives 



(*(r-,i)$(0,t)) 



r/2T 



(10) 



that is the correlations saturate (exponentially fast) for 
t > r/2. Note, however, that the precise details of the 
cross-over behavior for \t — r/2\ ~ tq are dependent on 
the form of F for < 77 < 1. 

Other CFT results. Many other general results may be 
found within the CFT formalism. For example, the 2- 
time correlation function 



($(r, *)$((),«)> 



e -irx(t+s) /4t 
e -TTxr/4T 
-7r|t-s|/4r 



where it is assumed that t, s, \t—s 



for r > t + s , 

for t — s < r < t + s , 

for?' < \t — s\ . 

(11) 

and r are all 3> tq . In 



the case where (<&) = 0, the first case gains an additional 
factor e -^b(t+s-r)/4To_ Note that the autocorrelation 

function with r = depends only on the time difference 
and does not exhibit aging in this regime. Another ex- 
ample is the f-point function in a semi-infinite chain, for 
which we find 



(*(r,t)) 



— TZXtj2To 
- ) — 7TXr/2TQ 



for t < r , 
for t > r . 



(12) 



The gaussian chain. The simplest exactly solvable model 
in which to study these effects is a general chain of cou- 
pled harmonic oscillators with a hamiltonian 



2 ^ 



2/2 



m 4>, 



•r+3 



(13) 



with the standard commutation relations imposed be- 
tween the <fi r and their canonical conjugates n r . The 2- 
point function can be found straightforwardly if tediously 
by integrating the Heisenberg equations of motion for the 
Fourier modes. The result is 



(14) 



cos(2ilkt)) 



BZ 



dk 



where fiofc and &k are the single-particle dispersion re- 
lations corresponding to Hq and H respectively, and the 
integral is over the first zone. The CFT results corre- 
spond to the case with ~ v\k\ as k — > 0, which gives 
(4> r (t)<f> (t)) ~ m (t - r/2v) for m (t - r/2v) > 1, while 
it vanishes for t < r/2v. Taking $ to be the conformal 
field e lq< ^ then reproduces all the CFT results before, with 
To ~ ttIq , x ex q 2 and F = 1. 

However, the advantage of this simple model is that 
the effects of other dispersion relations can be under- 
stood. In fact, the dominant contribution to (|14J) in the 
limits of large t and r comes from where, by a station- 
ary phase argument, the group velocity Vk = = r/2t. 
This suggests the picture given in the introduction of the 
correlations being due to left- and right-moving particles 
emitted at t = from closely separated points, moving 
with the group velocity. It also explains for the case of 
a semi-infinite chain why the relevant time scale is r/v 
rather than r/2v, since one of the particles is reflected 
from the end of the chain. For the energy density, pro- 
portional to (4> r +i — 0r) 2 , the asymptotic behavior is a 
constant, but the approach to this is dominated by the 
particles with the smallest group velocity, which, for a 
gapless lattice model, come from the zone boundary at 
|fc| = 7r. An explicit calculation^^ gives a correction 
~ £~ 3 / 2 cos(ri 7r t + 7r/4). Similarly, for a quench to a 
gapped lattice H, with fig > 0, the fastest particles cor- 
respond to a non-zero wave number. This gives rise to 
spatial oscillations in the correlation function. 
The Ising-XY chain. As an example consider the CFT 
describing the scaling limit of the Ising-XY chain with 
hamiltonian 

V-/ (1+7) » « ( l -l) „y„y hrr z\ ri c\ 
r=l 

where 7 is the anisotropy parameter and h is the external 
transverse field. It is well known that for any 7^0 
the model undergoes a phase transition at h — 1 that is 
in the universality class of the Ising model (defined by 
7=1). In that case F is known exactly [7j to be F = 
(s/l + rj 1 / 2 ± ^1-^/2) /y/2, where the upper and lower 
signs correspond respectively to fixed and free boundary 
conditions on the spins 0%. These correspond to quenches 
to the critical point h = 1 from ho < 1 and ho > 1 
respectively. Our general results should then apply, for 
example, to the correlations of the order parameter a x , 
with scaling dimensions x = i and Xb — \ for ho > 1 
and zero for ho < 1. 

The time-dependence of the transverse magnetisation 
(a*(t)) (which is in fact a piece of the energy density) 
has been studied in Ref. pj. For case of a quench from 
ho > 1 to h = 1 the results agree precisely with our gen- 
eral results of a constant asymptotic value approached 
by an oscillatory £~ 3 / 2 power law, which we have argued 
above is due to lattice effects. In Ref. |4fl the time aver- 



aged behavior of the spin-spin correlation function was 
calculated for quenches from ho = and infinity. For 
quenches to the critical theory with h = 1, the behavior 
was found to be purely exponential, in agreement with 
(HJ. In Ref. 0] the behavior of the correlation function 
in a finite system was studied numerically. While in this 
case there are complications from both the finite size and 
lattice effects, the main features are once again consistent 
with our general theory. 

We conclude with a discussion of the expected behav- 
ior in higher dimensions d > 1. For quenches from one 
disordered state to another outside the critical region 
the gaussian model results (|14f) (with J dk — ► J d d k) 
should apply. For t < r/2v m the 2-point function be- 
haves as at t = 0. At large times it decays exponentially 
with a scale set by the correlation length of H (although 
its Fourier transform differs in detail from the behav- 
ior in the ground state of H). The approach to this for 
t > r/2v m is however more complicated, since the wave 
fronts are no longer planar. A similar analysis may be 
applied to the dynamics of the Goldstone modes, or spin 
waves, in the quench to an ordered state. Above the up- 
per critical dimension (d = 3 if z = 1) mean field theory 
should also apply to quenches to the critical point. In 
d = 3, this gives ((/>(r,t)<p(0,t)} ~ r -i e -x(t-r/2v)/T for 
t < r/2v, and saturates at later times to t- independent 
value. Similarly we find, analytically continuing the re- 
sults in Ref. ^2] for a slab geometry, that for a critical 
quench from a disordered state, the energy density pa- 
rameter ((ft 2 ) decays exponentially. In principle it should 
also be possible to use the results in [13j for a slab geome- 
try in d — 3 — e dimensions, but this is very cumbersome. 
The calculation is simpler at large N 0] and once again 
shows exponential decay. However, a counter-example to 
this general rule appears to be given by the behavior of 
the magnetization in a critical quench from an ordered 
state. The mean-field profile in a slab geometry, given 
in involves elliptic functions, which, when continued 
to real time, give non-decaying oscillations with period 
~ To- We expect this result to be modified by the inclu- 
sion of fluctuations. 
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